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We describe a general method of calculating the fully differential cross section 
for the production of jets at next-to-leading order in a hadron collider. This method 
is based on a 'crossing' of next-to- leading order calculations with all partons in the 
final state. The method introduces universal crossing functions that allow a modu- 
lar approach to next-to-leading order calculations for any process with initial state 
partons. These techniques are applied to the production of jets in association with 
a vector boson including all decay correlations of the final state observables. 
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I. INTRODUCTION 



One of the most striking features of hadronic events is the appearance of 'jets' of hadrons. 
By use of a suitable experimental jet algorithm, the hadronic data may be organized into 
final states containing a definite number of jets. This defines the topological structure of the 
event for a given jet algorithm. Different jet algorithms or jet defining cuts can lead to a 
different number of jets being observed in a given event. As a result, jet cross sections depend 
on the procedure used to define an experimental jet. Nevertheless, at the experimental level, 
a jet is a perfectly well defined quantity since for a given jet algorithm each event contains 
a precisely determined number of jets. 

From the theoretical point of view, the jet algorithm plays an important role in selecting 
high momentum transfer events in which the soft radiation is removed by a cut on the 
minimal transverse energy of the jet. At the same time, individual hadron behavior is 
averaged out by replacing all hadrons within a cone of a given size by a single jet axis and 
jet energy. Because the hadronic information is averaged out, we can relate the hadronic 
jet axis and energy observed in the experiment to a jet axis and energy constructed from 
a parton shower calculated within perturbative QCD . The jet axis and energy obtained 
from the hadronic shower are thus modelled by the jet axis and energy obtained from the 
parton shower. This is a weak form of the parton-hadron duality theorem JTJ. Of course, 
non-perturbative hadronization effects are not predicted by perturbative QCD. Similarly, the 
contributions from the underlying event (at hadron colliders) are not included. 'Sensible' jet 
algorithms minimize these effects and allow a more direct comparison between theory and 
experiment. 

More precise theoretical predictions of jet cross sections are potentially valuable since new 
physics is often evident in events containing a specific number of jets. The corresponding 
QCD background is then the exclusive jet cross section containing the same number of jets. 
For example, the signal to background ratio for the top quark in the lepton + jets channel 
is improved by demanding that more jets be observed ||. 

The lowest order matrix elements for the two most prominent processes containing jets 
at hadron colliders, 

pp^rixjets, (1.1) 
pp W ± /Z + n 2 jets, (1.2) 

have been computed for n\ < 5 || and n 2 < 4 |4[] by making use of helicity amplitudes || , 
color decompositions |||7||1 and recursion relations to control the rapid increase in the 
number of Feynman diagrams as the number of partons involved grows. The jet cross section 
is then obtained by Monte Carlo integration of the phase space of the final state partons. 
This approach allows any experimental jet algorithm and acceptance to be applied, and one 
can study any distribution depending on the jet observables. It is important to note that 
at leading order, the jet is modelled by a single parton. The jet defining cuts are applied 
to this lone parton and the parton's direction and energy describe the jet's axis and energy 

Comparisons of lowest order QCD predictions of jet distributions with the data have 
proved reasonable, bearing in mind the fact that one is comparing a theoretical perturbative 
calculation with hadronic data. Generally, the shapes of infrared safe distributions are well 



predicted while the overall normalization is uncertain, due to a dependence on the unphysical 
renormalization and factorization scales hr and l fL2f . 

The addition of next-to-leading order effects produces three important improvements 
over a leading order calculation. First, the dependence on the unphysical scales and \ip 
is reduced so that the normalization is more certain. Second, we begin to reconstruct the 
parton shower. This means that two partons may combine to form a single jet. As a result, 
jet cross sections become sensitive to the details of the jet finding algorithm, particularly 
the way in which the hadrons are combined to form the jet axis and energy, and to the size 
of a jet cone. This sensitivity is also seen in experimental results. Third, the calculation 
becomes more sensitive to detector limitations, because radiation outside the detector is 
simulated. This can change leading order results considerably for quantities such as the 
missing transverse energy in events containing a W boson. 

From a more theoretical point of view, the presence of infrared logarithms in a generic 
perturbative QCD prediction implies that the expansion of physical quantities is not strictly 
a series in the coupling constant a s , but is rather a series in a s ln 2 y IR and a s \nyi R (as 
well as a s alone), where yj& is an experimental resolution (for example, the minimum jet 
invariant mass in e + e~ collisions). Thus even in the perturbative regime, the leading-order 
result — where the size of these logarithms is uncalculated — suffers from potentially large 
corrections which might spoil the applicability of perturbative QCD. In a next-to- leading 
order calculation, these logarithms are calculated explicitly, and thus one regains confidence 
in the applicability of the perturbative expansion. From this purist's point of view, a next-to- 
leading order calculation is necessary in order to understand whether a leading-order result 
is trustworthy. 

At next-to-leading order, the n-jet cross section receives contributions from virtual cor- 
rections to n-parton, and from real corrections in the form of (n + l)-parton, final states. 
Both contributions are divergent. The matrix elements for the virtual diagrams are infrared- 
divergent, while the real (n + l)-parton matrix elements are well defined. However, when 
the (n + l)-parton matrix elements are integrated over the allowed regions of phase space 
an infrared-divergent cross section is obtained. This comes about because the jet algorithm 
allows one of the partons to be soft, or for two partons within a jet cone to be collinear. 

In order to cancel these divergences explicitly, it is convenient to divide the (n + 1)- 
parton phase space into regions where (n + l)-partons are 'resolved' and regions where only 
n-partons are 'resolved' [IBJUI I5|,|IS| . For example, if the invariant mass of two partons, Sy, 
is smaller than some theoretical parton resolution parameter s m i n only one parton is resolved, 
while if S{j > s m ; n both partons are resolved. All of the divergences from the (n + l)-parton 
final state separate and are associated with the regions where only n-partons are resolved. 
These divergences can be cancelled directly against the virtual corrections to the n-parton 
cross section. With such a physical picture [[Kj it is straightforward to extend this method 
to deal with any number of partons in the final state. Indeed, using a color decomposition 
of the amplitude, one can write down a simple soft factorization for the sub-amplitudes, 
which in turn allows the construction of a universal set of functions summarizing the soft 
and collinear behavior of the matrix element of any colored particles, both before and after 
the cancellation of infrared divergences described above. Recently, we have described how 
this scheme may be applied to multijet cross sections in e + e~ annihilation This case is 
rather special since there can be no QCD radiation in the initial state and all divergences 



therefore reside in the final state. In the paper, we wish to extend this method to processes 
which involve partons in the initial state. 

In Sec. U we show how to set up the next-to-leading order calculation of the matrix 
elements for e + e~ — *> n partons and how they may be 'crossed' to give the cross section for 
pp — > V + (n — 2) partons where V = W ± , 7*, Z. In order to do this it is convenient to 
introduce (in addition to the universal soft and collinear functions mentioned above) a set 
of universal crossing functions which multiply the lowest order cross section. These crossing 
functions are essentially convolutions of the input structure functions with the Altarelli- 
Parisi splitting functions. Section [III] deals with the derivation of the crossing functions 



which enable us to cross final state partons to the initial state. In Sec. [TV], we construct 
explicit Monte Carlo programs for, 

pp -> V + 0, 1 jets -> ££' + 0, 1 jets, (1.3) 

at next-to- leading order. The jet algorithm may then be applied directly to the n and (n + 1) 
parton final states. All dependence on the unphysical parton resolution parameter s min 
cancels numerically . The cross section is fully differential in all jet and lepton observables 
and therefore differs from calculations of the W transverse momentum distribution at 0(af) 
p2|p3| or of the single jet inclusive transverse momentum distribution at 0(a? s ) |T3] , |24|] . 
Equivalent techniques have been applied to pp — > W ± + jets |2S|] and pp — > 2 jets [T2(J at 
next-to-leading order. Finally, we summarize our results in Sec. |V]. 



II. CALCULATIONAL ORGANIZATION 

In this section we give a general overview of the manner in which the calculation is set up, 
without going into details, which we shall present in subsequent sections. The organization 



minimizes the computational effort while retaining the standard MS prescription ||27|| . For 
example, the cumbersome (i-dimensional squaring of the matrix elements is avoided. Fur- 
thermore, the introduction of universal crossing functions will allow us to obtain the initial 
state parton cross sections from the all-outgoing cross section. These techniques depend 
crucially on the universality of the QCD soft and collinear radiation patterns. 

An efficient way of organizing next-to-leading order calculations for all-outgoing parton 
processes was given in ref. [ITH]. The basis of this method is the use of ordered amplitudes 



associated with each color structure rather than the full matrix element. For example, let 
us consider the decay of a vector boson into a quark-antiquark pair with n accompanying 
gluons. The full squared amplitude is obtained by summing the squared ordered amplitudes 
over all permutations of the gluons 



\M v (V->qq + ng)\ 2 ~ £ | A(q; g u . . . , g n , q) | 2 . (2.1) 

Perm. 



For simplicity, we keep only the terms at leading order in the number of colors. See Sec. IV for 



a full discussion of the subleading terms. Of crucial importance is the fact that the ordered 



amplitudes exhibit QED-like factorization [23] in the soft and collinear limits |[21||. This 



forms the basis of our method and allows the integration over the singular (or unresolved) 
parts of phase space without calculating the hard matrix element explicitly, 



d Psofa I collinear 



A(Q'i9i,-- -,9^+11^) R(q; 9i, ■ ■ ■ , 9n, q\Smm) \ A(q; g u . . . , g n ; q)\ (2.2) 
where the soft and collinear regions of phase space are defined by the invariant mass cut s m i n . 



At this point everything is done in an arbitrary number of dimensions fl29| , |30|j . However, we 
never have to calculate the hard amplitude squared explicitly since this result is obtained 
without any detailed knowledge of the hard process. The next step is to calculate the virtual 
corrections to the squared matrix elements which have the generic form, 



\M V (V -»• qq + n g)\ 2 y = ^ V(q; g x , . . . , g n ; q) \A(q\ g u . . . , g n ; q)\ 2 + T{q\ gi, . . . , g n ; q), 



2 E 1 

Perm. 

(2.3) 



where V(q; gi, . . . , g n \ q) is the singular part proportional to the tree level ordered amplitude 
and T{q\ g\, . . . , g n ; q) is the remaining finite contribution. This can be immediately com- 
bined with the unresolved phase space contribution Eq. [2.2| to give the finite next-to-leading 
order squared matrix elements, 

\M v (V^qq + ng)\ 2 F ~ £ ([1 + K] \A(q; g u . . . , g n] q)\ 2 + F(q;g u . . . , g n] qj) . (2.4) 

Perm. 



Note that due to the Kinoshita-Lee-Nauenberg PT],P2[ theorem the combination of the phase 



space factor and the virtual factor, K — R + V, is finite. As a direct consequence we can 
now perform the squaring and summation over the polarizations in 4 dimensions using the 
standard techniques developed for evaluating complicated tree level amplitudes (helicities 
and recurrence relations |§). While V can be calculated in a process- independent manner 
jnj| , the finite remainder of the virtual correction T needs to be calculated on a process-by- 
process basis. The general structure is process-independent, and it is in this sense that the 
fC factor is universal. 

In order to generalize the framework above, used in e + e~ collisions, to hadronic collisions, 
we must include initial state partons in the calculation. One useful property of lowest order 
matrix elements is that of "crossing" . In other words, the matrix elements for V — > qq + ng 
are related to those for the crossed processes, 

qq — ► V + ng, 

qg^V + q + (n-l)g, (2.5) 

qg -> V + q + (n - 1) g, 

99 ~>V + qq + (n-2)g, 

by reversal of the momentum and helicity of the crossed particles. The fully differential cross 
section at leading order in the collision of hadrons Hi and H 2 , 

H x + H 2 -^V + n partons, (2.6) 



is. 



d<r Hl H 2 =y2fa 1 { x i)fb H2 ( x 2)da^(x 1 ,x 2 ) dxidx 2 . (2.7) 



ah 



Here fa{x) is the probability density of finding parton a in hadron H with momentum 
fraction x and, 



do^ b °(xi, x 2 ) = ab \M-ab^ dP(ab^V + npartons), (2.8) 

" $ab 

where $ a fc is the appropriate spin and color averaging factor and dP the V + n parton 
phase space where all parton pairs satisfy s i? - > s min . The matrix elements for ab — > V + n 
partons are denoted | and are related by crossing to \Mv\ , Eq. gj. 

We have already discussed how the next-to-leading order matrix elements for e + e _ — > 
n partons can be written in an explicitly finite way using the parton resolution parameter 
•Smin- We now wish to extend this to processes involving partons in the initial state while 
maintaining the crossing properties of lowest order. In order to achieve this, the next-to- 
leading order hadronic cross section must be defined by, 

da Hl H 2 = Y,^a 1 (xi)J r b H2 {x2)dcr^ b LO {x 1 ,x 2 ) d Xl dx 2 , (2.9) 

ab 

where is the "effective" next-to-leading order structure function and do^ b ° the "crossed" 
analogue of the finite next-to-leading order partonic cross section. This cross section can be 
expanded as a series in the coupling constant, 

d*l LO = dog + a s d5a^° + 0{a% (2.10) 

where we have extracted the coupling constant from the finite crossed matrix elements. Note 
that a s is evaluated at the renormalization scale fiR. Similarly, after mass factorization, the 
effective structure function may be written as, 

Fa(x) = fa(x^F) + a s C?(x^ F ) + <D(a 2 s ), (2.11) 

where hf is the factorization scale. Both f^(x,fip) and the crossing function C^^x^hf) are 
finite. Once again a s is evaluated at the renormalization scale. In principle one could evaluate 
a s at the mass factorization scale, however, provided a s log(/i|.//i^) <C 1, the difference is of 
0(a 2 s ) and can be ignored. For a detailed derivation of the structure of the crossing function 
we refer the reader to Sec. llTT. 



Inserting these definitions back into Eq. fZM and expanding up to 0(a s ) we find, 



do Hl n 2 — ^2 



ab 



/f^O/f'fo) {dag{x 1 ,x 2 ) + a s d5o^°{x 1 ,x 2 )} 



+ a s {C^( Xl )f b H \x 2 ) + f^{ Xl )C^{x 2 )} da L ab °(x u x 2 ) + 0{ 



dx\dx 2 . 

(2.12) 



For simplicity, we have suppressed the dependence on the renormalization scale in the cou- 
pling constant and the factorization scale in the structure and crossing functions. 

The crossing function receives two contributions which both stem from the fact that we 
consider two partons to be unresolved when their invariant mass is smaller than the parton 



resolution parameter s min . Firstly, we cannot distinguish between a single initial state parton 
and a parton which emits collinear radiation such that the invariant mass of the collinear pair 
is smaller than s mm . This implies that part of the initial state collinear radiation is removed 
from the hard scattering and absorbed into the effective structure function. Clearly, this 
contribution depends on s min and therefore so do the crossing functions (and also T^Y This 
term is a convolution of f£* with the Altarelli-Parisi splitting function, P a _, c (z) ||33|| . 

The second contribution arises from crossing a pair of collinear partons with an invariant 
mass smaller than s min from the final state to the initial state. In principle we should remove 
this contribution from dSa^°, however in order to preserve the structure of Eq. |2i], we 
subtract this contribution from the parton density function. This is possible because we 
cannot distinguish the two parton incoming state with invariant mass smaller than s m i n 
from a single incoming parton. 

Both of these contributions are divergent and schematically, 



x Z \ Z , 



dzP a ^ c {z) 



(2.13) 



A more precise formulation of the crossing function is given in the next section including all 
c/-dimensional factors. After mass factorization the crossing function is rendered finite, 



a 



//.scheme ( 



X,fl F ) 



N 
2^ 



AZ(x) log 



/4 



+ B, 



//.scheme , 



X) 



(2.14) 



Although is scheme independent, does depend on the mass factorization scheme and 
therefore so does C„ . Explicit forms for these functions in the MS scheme are given in the 
next section. 

The overall cross section cannot depend on the unphysical parameter s mm . When the 
contribution from H\ + H 2 — > V + (n + 1) partons where all partons are resolved is included, 
the s m in in the logarithm is replaced by an energy scale defined by the experimental cuts. In 
a numerical computation, one would also force the factorization scale to be determined by 
the experimental cuts; the argument of the logarithm would then be of 0(1) and the contri- 
bution from A^ would be small. If this were not true, the logarithm would be large so that 
a s A^(x) log (-^exp/ Pf) — fa( x )i perturbation theory would break down and a resummation 
of the leading logarithms would be necessary. 



III. DERIVATION OF THE CROSSING FUNCTIONS 

In this section we derive explicit formula? for the crossing functions C^(x,^lf) as de- 
fined in eqs. [2.9^2. 12j . First, we derive the initial state collinear phase space behavior in 



a parametrization suitable for our parton resolution parameter. We then reformulate the 
standard collinear matrix element factorization in the ordered amplitude language. Com- 
bining these results enables us to derive the universal crossing functions which after mass 
factorization yield finite crossing functions. 



A. The initial-state collinear behavior of phase space 

First consider the production of a heavy object, Q, by the collision of two massless 
particles with momenta p a and p h . The <i-dimensional phase space measure, including the 
flux factor, is given by, 

1 2n 
-dP\a + h^Q) = —5(s ah - Q 2 ), (3.1) 



2s a ^ 2s a /j 

where s a h = (p a + Ph) 2 - This extends straightforwardly to the production of any number of 
particles (massless or otherwise) with momenta p±, . . . ,p n by use of the relation, 

dO 2 

dP d (a + /i->l + -- - + n) = dP d {a + Q)-^-dP d (Q -> 1 + • • • + n). (3.2) 

2tt 

Next consider the phase space for the production of a massless particle with momentum 
p u in association with Q from the collision of two massless particles with momenta p a and 
P P , 



1 7i I Ti I r | || | -| — (2— d/2) 

dP d (a + p^u + Q) = (27r) 2 - d lSaul lSpul 



?au | | J pu | 



Sap 



<ti}d--3 r, _ | | _ | | 

g 0\Sap | Sau| \S p u\ 



(3.3) 

where we integrate over the invariant mass of Q and the polar angle with respect to p p by 
using |s au | and \s pu \, as well as integrating over the (d — 3) azimuthal angles relative to the 
direction of p p . 

The region where momentum p u is collinear with momentum p p is defined by, 

In this region we introduce the hard momentum p h which is the amount of the parent 
momentum p p remaining after the emission of the unobserved collinear momentum p u such 
that, 

Ph %Ppi $ah ZS ap i 

p u = (1 - Z)p p , \s au \ = (1 - Z)s ap - (3.5) 

In this limit the phase space factorizes, 

' dP d (a + p->u + Q) -> dP d oX (p ^u + h)x ^—dP d {a + h^Q), (3.6) 



2s ap 2s ah 

where, taking d = 4 — 2e, we find, 

dP^ 2e (p ^u + h) = ^^—-zdzd\ Spu \ [(1 - z)\s pu \r ■ (3.7) 

The square bracket contains the necessary factors to regulate the poles in the matrix elements 
in (1 — z) and s pu . 

Combining these results we find that in the collinear limit, the full phase space measure 
of interest factorizes as follows, 

dP d {a + p^u + 2 + --- + n) = dP d ol {p -> u + h) x —dP d (a + h^2 + --- + n). (3.8) 

Sah 



B. Behavior of matrix elements 



The matrix elements also undergo a collinear factorization when one of the final state 
partons is collinear with one of the initial state partons. Take the case where an initial state 
parton p splits into partons u and h (which participates in the hard scattering) as in Eq. 
then, for each ordered amplitude, 



A(...,p,u,n,...) 



*uh 



A(...,h,n,...) 



where, 



*uh 



'g 2 N\ 1 P hu ^ p {z) 



(3.9) 



(3.10) 



J pu I 



Note that the quantum numbers of the unobserved parton u are determined by the quantum 
numbers of the parent parton p and the hard-process parton h. 

This is very similar to the factorization that occurs when two final state particles are 
collinear. In this case, when parton a (which participates in the hard scattering) splits into 
a final state collinear parton pair 1 and 2 then, 



A(. 



,1,2, 



Cp 



A(. 



where, 



(3.11) 



Sl2 



(3.12) 



and, 



Pi = Zp a 



P2 



Z)p a - 



(3.13) 



The different averaging factors for initial- and final-state quarks and gluons have been taken 
into account, however, we do not sum here over different flavors of quarks participating in 
the hard process. 

As before, the splitting functions may be either in the conventional scheme (all particles 
in d-dimensions) ]33|] or in the 't Hooft Veltman scheme (only unobserved particles in d- 
dimensions) ||29|| . In the conventional scheme, the splitting functions are given by, 



= ( 1 - ^ ) p< 



1-z 

■ig->q( z ) : 



: 1 — z .„ 
+ + z(l 

z 



l + z 2 



:(1-Z) 



1 - Z 



gg~^g 



» 



jy-Pqq-*g( z ) 



2 

iV 



z 2 + (l- z) 
1-e 



(3.14) 



For the splitting functions in the 't Hooft Veltman scheme, see Ref. P9p . 

One difference from pure final state singularities is that the initial state parton is always 
hard — there is always a minimum value for z imposed by demanding that a hard scattering 



takes place. On the other hand, the upper bound on z is still determined by the requirement 



that parton u is collinear but not soft. In other words, s l 



> s r 



where n is the neighbouring 



hard parton in the ordered amplitude (see eq. |3.9|) . 

For the g — > gg process, there will be contributions from two ordered amplitudes, 

2 2 

A(. . . , m,p, u, n, . . .) + A(. . . , m, u,p, n, . . .) — > A{. . . ,m,h,n, . . .) 



(3.15) 



where the order of the other hard partons in the ordered amplitude is preserved. The upper 
limit on z will be different in each case since the requirement that gluon u be unobserved 
depends on the adjacent momenta. Note that in the final state case, each ordering counts 
equally, however, the Bose symmetry factor takes this into account. For processes involving 
collinear quarks (antiquarks), only one ordering will contribute. Note that only ordered 
amplitudes where p and u are adjacent contribute in the collinear limit. If they are not 
adjacent as in Eq. |3.15 the collinear limit gives a contribution of the order of the parton 



resolution cut s min , which is therefore negligible. This property is very useful since it avoids 
overlapping divergences for a given ordered amplitude and this makes partial fractioning to 
isolate the divergences unnecessary. 



C. Behavior of the cross section 

In this subsection we derive exact expressions for the crossing functions C^ix) 
alently the effective structure function T^[x) as defined in Eqs. |2.9| - |2.12| . We will consider 
the generic process of scattering of partons a and h to form an arbitrary final state with an 
invariant mass (e.g. partons only, vector boson plus partons, etc.). The leading order 
cross section for the production of a vector boson plus partons is given in Eq. |2?7]. Cross 
sections for other final states are given by similar formulae. The next-to-leading order cross 
section is defined in Eq. |2.9| , or in its expanded form in Eq. |2.12 . 



The first step in the derivation of the crossing function is to consider the initial-state 
collinear radiation contribution to the next-to-leading order cross section. Consider the 
splitting of a parent parton p to a (unobserved) collinear parton u, and a parton h partic- 
ipating in the hard scattering: p — > uh, where the invariant mass \s pu \ < s m j n so that this 
configuration is indistinguishable from the leading order configuration where parton h comes 
directly from the hadron. This contributes to the next-to-leading order cross section and 
using Eqs. |3]7| and we find 

^initial = fa 1 ^!) {fp'"^^^^^ + h)5(x 2 - zy)dy} da^(x u x 2 )dx 1 dx 2 , 

ahp 

(3.16) 

where, by definition, the momentum fraction x 2 carried by parton h is given by the momen- 
tum fraction y of the original parton p multiplied by the energy fraction remaining after 
radiating the unobserved parton u. There is an implicit integration over z contained in 
the collinear phase space factor, see Eq. [3.7| . Comparing Eq. |3.16| with Eq. |2.12| gives the 



contribution of the initial state radiation to the crossing function, 

«^ 2 initiai(* 2 ) = E fp\y)cV uh dPUp - u + h)6(x 2 - zy)dy. (3.17) 



Using Eq. [3.10| and the collinear phase space factor of Eq. |3]7| gives, 

<*> - (£) ^ o' i £ \ r i* 1 - -)-A^M? (? 

iV\ 1 /47r/i 2 \ e 1 ^ f 1 dz H fx 2 



E/ -//(^U^M, (3.18) 



27T/ T(l - e) \ s min y e p Jx2 z 

where we have integrated y over the delta function and s pu using the constraint |s pu | < s min . 
Here \i is an arbitrary scale introduced to keep the strong coupling constant a s = g 2 fi~ 2t / 4tt 
dimensionless in ^-dimensions. The upper boundary on the z integral is determined by the 
constraint that the unobserved parton u is not soft with respect to its neighbouring parton 
n. In other words, \s un \ = (1 — z^s^l > s min . Explicitly this gives 



z< l__Vm l_^ 2 . (31g) 



Note that Shn is only defined because of the use of the ordered amplitudes and is different 
for each ordering. 

Looking at the definition of the splitting functions, Eq. |3.14| , we see that J q ^ g and J g -> q 



do not depend on the upper boundary on z (up to negligible corrections of 0(s min )). This 
is due to the absence of a singularity in the limit that the quark or anti-quark becomes 
soft. In contrast, the J g ^ g and J q ^ q functions do contain a soft singularity, arising from 
the limit where the gluon becomes soft, and therefore do depend on z?,- In order to write 



these contributions to the crossing functions in the second form in Eq. |3.18| , we use the ( )+ 
prescription defined by, 

(F(z)) + = Km Ul - z - P)F(z) - 5(1 - z - /?) F(y) dy\ , (3.20) 



such that, 
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provided that g(z) is a function well behaved at z — 1. 



The functions J p -+h are thus given by, 
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(3.25) 



The next step is to correct for the fact that we have crossed a final state collinear cluster 
to the initial state. As explained in Sec. |T| this is done by subtracting the collinear factor 
resulting from the splitting h — > up integrated over the final state collinear phase space (see 
TE for a detailed derivation). The contribution to the next-to leading order cross section is 
given by, 



ahp 



d 

col, final 



[h — > u + p)\ da^h (x\, x 2 )dxidx 2} (3.26) 



giving, 
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Note that the parton density function is associated with parton h rather than with parton p. 
The integration boundaries of the z integral are again defined through the requirement that 
the hard partons are resolved. For each ordered amplitude z\ and z 2 are given by demanding 
that the invariant mass of parton u with both its neighbors in the particular ordering is 
larger than the parton resolution cut s m i n so that u is not soft. In the conventional scheme, 
the final state integrals over the splitting functions, I pu ^h, are given by, 
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(3.29) 



We can now define the crossing function as a convolution integral involving the parton 
density function and a crossing kernel X p ^h{z) which is obtained by subtracting the final 
state contribution given by Eq. |3.29| from the initial state contributions of Eq. |3.25 
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where the crossing kernel for specific processes is given by 



iV 



1 



'47T/! 2 



2tt/ T(l - e) \ s min / 
x ( J g ^ ff (z, zi) + J 9 ^ fl (z,z 2 ) 



z 1 ,z 2 ) + ^-I qq ^ g (0,0) 



6(1 -Z) 



N 



( 4:71/1 



2\ e 



X 



2W r(l -e) V s 

( z 



\l-e 



+ 



J(l-^) 1+ 1 + 

UN- 2n f 



6iV 



/V 67 5n/\\ 



= -(i)f(^) O' I M- 2 *> - Wo. *W - .)) 



l 



s 

47T/i 



1 

X - 

e 



2tt/ r(l -e) 
l( l + z 2 



mm 

1 



AT2 



7T 



2 V[(l-^) 1+e 



t(1 " S) + I- f 6-4 S(l - Z 



We see that the dependence on the boundaries exactly cancels, making the crossing function 
independent of the hard process. The other two functions, X g ^ q (z) and X q ^ g (z), do not 
receive contributions from the final state crossing, 
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(3.31) 



Here, we have replaced endpoints that do not contribute with a zero (that is, we have simply 
dropped contributions of 0(s m in))- 

We have now derived the process independent crossing function. They still contain the 
mass singularity which has to be removed by the mass factorization prescription. In fact in 
the language we have developed here the mass factorization is done very easily as is shown 
in the next subsection. 



D. Mass factorization 



The only physical, and therefore finite, quantity associated with resolved partons is the 
effective structure function J-'^x) as defined in Eq. [2 . 1 1| . Conventionally, the parton density 
function is made finite by renormalizing the parton density function at the factorization scale 
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This is very similar to coupling constant renormalization. The 0(a s ) term is subsequently 
absorbed in the crossing function 



tf(x) = tf(x) + a.CZ(x) 
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with 



(3.34) 



where the mass factorization counter function absorbs the divergences in the crossing func- 
tions. Note that the effective structure function, ' , is left unchanged by the mass fac- 
torization and is in fact independent of the factorisation scale. However for a cross section 
calculated at fixed order in perturbation theory, we have to expand the effective structure 
functions explicitly and neglect terms of 0(a 2 ) as was done in Eq. |2.12| . This makes the 
fixed order cross section factorization scale dependent since we have to neglect the term 
^( MF )C^(x 1)Mj ,)Cf 2 (a; 2) /x F ). 

The mass factorization counter functions R p ->h at the factorization scale \x F are given by, 
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(3. 

is the d — 4 part. The function fl c }^ ae (z) is the scheme dependent mass factorization term 



where the four dimensional part of the splitting function is given by P^ b ^ c (z) and P t 

xschei 



chosen such that f¥^ h (z) = 0. The strong coupling constant in (3.32) is evaluated at the scale 
/i, which through coupling constant renormalisation is identified with the renormalisation 
scale. Other choices of the scale are possible, however, provided a s log( fi 2 R / fi F ) <C 1, the 
difference is of 0(a 2 ) and can be ignored. Indeed, this condition is necessary in order to 
prevent the appearance of large logarithms. In practice, [ir and hf will usually be chosen 
equal, but if they are not, the ratio /ir/ '/if should be small. 

Combining the unrenormalized crossing functions of Eqs. ??- 3.31 with the counter func- 
tions of Eq. p.35| gives us the finite, renormalized crossing functions in the MS scheme, 
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where the arbitrary scale [i has canceled and, 

A% (x, n F ) = Af^ h (x, fi F ), 



(3.36) 
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The finite scheme independent functions A^ h (x, hf) are given by, 
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and the scheme dependent functions -B5^ s (x) by. 
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With these formulae we can calculate the crossing functions for each set of given parton 
density functions. These crossing functions are independent of the hard process. It is now 
straightforward to evaluate Ah(x, fxp), Bh(x,^ip) and C/^x, /xf) numerically for a given set 
of input parton density functions in a particular scheme. We use the MS scheme, and in 
order to give some idea of the size and relative importance of the crossing functions we use 
set Bl of Ref. as input proton MS parton density functions. Furthermore, we focus on 
the crossing functions associated with valence up quarks and gluons. The distributions for 
down valence quarks show a similar behaviour to the up valence quarks while the sea quarks 
are related to the gluonic distributions. 

First, we show the x dependence of the crossing functions for valence up quarks at a fixed 
scale up = 25 GeV in Fig. 1. In order to illustrate the s m j n dependence of C u , three values 
of s min have been chosen, s min = 1, 10 and 100 GeV 2 . The first two values are typical of 
the s min chosen in practical applications (see Sec. [TVp and are values where the systematic 
uncertainty in evaluating the cross section is of the same order as the uncertainty introduced 
by approximating the matrix elements at small ,s mm . Although C u does explicitly depend on 
this unphysical parameter, this dependence is balanced by a growth of the bremstrahlung 
contribution to the next-to-leading cross section. Once this cancellation takes place, s m i n is 
replaced by a scale of order of the experimental cuts. For jets with Et > ETmin ~ 0(15 GeV) 
and a jet-jet separation of AR ~ 0.7, this scale is of order E^ min AR 2 ~ (9(100 GeV 2 ). The 
curves with s min =100 GeV 2 are representative of such a scale and give some indication of 
the contribution to the physical cross section. As shown in Fig. 1(a), both A u and C u are 
negative for some values of x. Although this seems somewhat strange, A u and C u are not 
directly interpretable as physical distribution; only the complete next-to-leading order cross 
section as defined in Eq. |2.12| is expected to be positive (so long as higher-order corrections 
are not too large). It is worth noting that although A u and B^f s are roughly similar in size, 
the contribution of A u to C u is enhanced by log(s mm //^) which can be large. Therefore, 
particularly for small s mm , the shape of C u is dictated by the scheme independent function 
A u , while B u is only important as s m j n — > fip. As mentioned earlier, a more physical quantity 
is the effective structure function T u defined in Eq. |2.11| . This is shown in Fig. 1(b) for the 
same three values of s min as is the ordinary parton density function f u . We see that at 
large x, T u is enhanced relative to f u , while at small x there is a depletion. Furthermore, 
as s m ; n — > fi F , T u approaches f u . Note that for very small factorization scales such that 
ji 2 F < s m i n , then T u is depleted at large x and enhanced at small x. 

Fig. 2 shows the x distribution for the gluonic crossing functions. These crossing functions 
receive contributions from both g — > gg and q — > gq splitting functions, and, due to the soft 
gluon poles that are present, A g and B g both grow at small x. As a consequence, C g is 
negative in this region. However, for large s mm , there is a significant cancellation between 
A g and B g so that C g is less singular. This is reflected in Fig. 2(b) where T g and f g are 
shown for the gluon. At small x there is a dramatic softening of the growth of the gluon 
density function. This depletion is entirely consistent with the depletion of the up valence 
distribution discussed above. Similarly, at large x, there is a small enhancement. 

It is also instructive to study the scale dependence of the crossing functions. This depen- 
dence is present in the input parton density functions and hence A and B and through the 
log(s min //i 2 7 ) factor multiplying A in Eq. |3.36| . T contains an additional \ip dependence from 
the strong coupling constant evaluated at the factorisation scale which we take to be equal 



to the renormalisation scale. Fig. 3 shows the [if dependence for the up valence density 
functions at x — 0.05. A u is a slowly decreasing function, much smaller than B u while both 
A u and B u are positive for /ip < 1000 GeV. As a consequence, for s min < /i F , C u is relatively 
small due to a cancellation between the two terms. This is not the case for smaller scales 
where C u can be quite large. Fig. 3(b) again shows the ordinary parton density function f u 
with the effective structure function T u . For large fip, where C u is small, f u and T u are very 
similar in size, while at small scales we can see very large differences. 

The scale dependence of the gluonic crossing functions is shown in Fig. 4. Unlike the up 
valence quark case, A g becomes negative at a scale of fip — 7 GeV. As a consequence, A g 
and the rather large B g combine together coherently to form C g which grows logarithmically 
at large scales. Furthermore, the effective structure function T g is always significantly larger 
than the ordinary parton density function f g . 



FIGURES 



FIG. 1. The valence quark density functions (a) A u , B u and C u and (b) T u and f u as a function 
of x for \xp = 25 GeV. C u and T u are shown for s m j n — 1, 10 and 100 GeV 2 . 



FIG. 2. The gluon density functions (a) A g , B g and C g and (b) T g and f g as a function of x 
for up = 25 GeV. C g and JF S are shown for s m i n = 1, 10 and 100 GeV 2 . 



FIG. 3. The valence quark density functions (a) A u , B u and C u and (b) T u and f u as a 
function of the factorisation scale fip in GeV for x = 0.05. C u and T u are shown for s m ; n — 1, 10 
and 100 GeV 2 . 



FIG. 4. The gluon density functions (a) A g , B g and C g and (b) J- g and f g as a function of the 
factorisation scale \ip in GeV for x = 0.05. C g and T g are shown for s min = 1, 10 and 100 GeV 2 . 



IV. APPLICATIONS 



We turn next to the construction of an example of a next-to-leading order cross sections 
for jet production at hadron colliders. Let us focus on the process, 



V — > qq + n g, 
for which the lowest order matrix element is given by, 

M v (Qu 1, • • • , n; Q 2 \P) = 5„(Qi; l,...,n; Q 2 )V». 



(4.1; 



(4.2) 



Here V 1 represents either the vector boson polarization vector, the lepton current which 
created the vector boson, or the leptonic decay products of the vector boson, while is the 
hadron current. Both currents depend on the particle helicities which we suppress throughout 
and the particle momenta which we denote by P for the vector boson, Qi, Q 2 for the quarks 
and Ki, . . . , K n for the gluons. These momenta satisfy the momentum conservation relation, 



P M = Q? + Q 2 + K{ + ■ • ■ + K£. 



(4.3) 

, a n and the quarks 



In addition, the hadron current depends on the colors of the gluons a±, 

Cl,C2- ^ 

The hadron current may be decomposed according to the different allowed color 
structures H0,|J, 



SJQx;l,...,n;Q 2 



ieg n £ (T^---T a -) ClC2 S,(Q i; l,...,n;Q 2 ), 

P(l,...,n) 



(4.4) 



where represents the colorless ordered amplitude in which the gluons are emitted in an 
ordered way from the quark line. Note that the prefactor associated with each is also 
ordered according to the color of the gluons. These color factors form a complete basis and 
therefore each is gauge invariant with respect to the gluons. Although the full hadron 
current is invariant under permutations of the gluons, the ordered amplitude is not. This 
property is recovered by summing the ordered amplitudes over the n! gluon permutations, 
P(l,...,n). 

One advantage of using this color decomposition is that the squared matrix elements 
summed over helicities and colors have a very systematic structure, 
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where n > 1 counts the number of gluons. On the right hand side, we have expanded in 
the number of colors. The terms subleading in the number of colors are related to matrix 
elements with abelian couplings. For example, when n — 2, 
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where, 



S M (Qi; 1, 2; Q a ) = «5 M (Q i; 1, 2; Q 2 ) + 5 M (Q i; 2, 1; g 



2^ 



(4.7) 



and the contribution from the triple gluon vertex drops out. The Feynman graphs contribut- 
ing to S^(Qi, 1, 2; Q 2 ) are therefore those for V — > gg + 2 7 or — > gg + (7 + 7. 

The ordered amplitudes also have special properties in the soft gluon or collinear parton 
limits which allow us to isolate the singular regions using the parton resolution parameter s m i n 
^T[|15|]. These divergences are proportional to the lowest order squared ordered amplitudes 
as are the virtual divergences. Therefore we can combine them directly and, due to the 
Kinoshita-Lee-Nauenberg |f31| , |32"| theorems, obtain a finite result after the usual coupling 



constant renormalization. The finite next-to-leading order matrix elements can be written, 
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(4.8) 



where J-{Qi] 1, . . . , n; Q 2 ) is the finite contribution from the virtual graphs. The dynamical 
ordered K factor is given by, 



K.(Q 1 ;l,...,n;Q 2 ) 




+ a s ([i R )b n\o. 



(4.9) 



where the sum runs over all (n + 1) color-connected pairs, that is, ij = Q\l, 12, . . . , nQ 2 , and 

where b is the one-loop coefficient of the QCD beta function. Note that fC(Qi, l,...,n; Q 2 ) 

depends explicitly on the parton resolution parameter s m i n . The renormalisation scale is 

the scale at which the MS counter term is subtracted. For vector boson decay all Sy > so 

that Q(sij) = 1. However, when we cross partons into the initial state this will no longer 

be true and it is necessary to maintain the explicit analytic continuations of the log 2 terms. 

Up to this point, we have continued both matrix elements and phase space into d = 4 — 2e 

using dimensional regularization. In Eq. |4.9| , we now see that all singularities have cancelled 

explicitly and we may therefore take the 4-dimensional limit. This means that the squared 

2 

ordered amplitudes, S^V* , may be evaluated in 4-dimensions and not in <i-dimensions, 
thus simplifying the calculation dramatically. Similarly, it is not necessary to extend the jet 
algorithm to (i-dimensions as in the work of Ellis, Kunszt, and Soper |24 . 

the 



Keeping all orders in the number of colors presents no problems. For example, 
effective matrix elements for V — ► gg + 2 g at next-to-leading order are given by, 
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(4.10) 



where K,{Qi\ Q 2 ), JC(Qi, 1; g2) an d JC(Qi, 1, 2; g 2 ) are given by Eq. £Rj with n = 0, 1 and 2 
respectively. 



A. ffiiJ 2 -» F + jets 



As a first application of p. 12| and the next-to-leading order crossing approach, let us 
consider the production of a vector boson in hadron-hadron collisions followed by the decay 
of the vector boson in the absence of jets. The relevant parton-level processes are, 



qq -> V, 

along with the bremstrahlung processes, 

qq -»• V + g, qg->V + q 
for which the generic cross section is given by, 



V + q, 



da ab (xi,x 2 ) 



2s afe 



\M ab \ dP(ab^V + 0, lpartons) 



The spin and colour averaging factors, $ a f, are given by, 
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(4.12) 
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The lowest order matrix elements for these processes are related to those for Eq. [O] with 
n = and 1 through the usual crossing relations. In other words, the momenta and helicity 
of crossed particles are reversed. For example, 



M qg (Qr, 1, . . . , n; Q 2 \P) = M v (-Qi; 1, • • • , n; -Q 2 \ - P) 



and, 



M gg (Qi;l,...,n;Q 2 \P) 



M v {Qi;-l,...,n; -Q 3 \-P) 



(4.15) 



(4.16) 



An explicit form for these matrix elements and hence |«S M (Qi; Q2)V >J '\ 2 and {S^Qi; 1; Q 2 )y^\ 2 
using spinor language is given in Appendix A of |l^] . Full details of how crossing affects the 
spinors is given in Appendix E of Ref. j35] . 



At next-to-leading order, the finite effective matrix elements for Eq. fO] with n = 0, are 
given by, 



M v {Q l] Q 2 \P)^ F = e 2 N (l-i 



(4.17) 

where, because of our choice for the assignment of the finite pieces between /C and 

^(Q i; Q 2 ) = 0. (4.18) 

The dynamical /C factor is given by Eq. [D| with n — 0. Because of Eq. |2.12| , we may cross 
this in exactly the same manner as the tree level matrix elements of Eq. |4.15| . The only 



subtlety is in the analytic continuation of the log 2 terms in /C which we have written out 
explicitly in Eq. |4.9| . 

It is now straightforward to construct a Monte Carlo program to evaluate the fully dif- 
ferential cross section numerically. In particular, the vector boson decays are easily included 
which allows experimental cuts to be placed directly on the observed leptons. It is important 
to note that the phase space is restricted to regions where all partons are resolved. In other 
words, any pair of partons must have an invariant mass larger than the parton resolution 
parameter, |sy| > s m ; n . What this means in practical terms is that the bremstrahlung contri- 
bution to the cross section grows as log 2 (s min ). This is balanced by the explicit — log 2 (s min ) 
in /C such that the total cross section should be independent of the unphysical s min pro- 
vided (a) that s m ; n is small enough that the soft and collinear approximations are valid and 
(b) that s m in is not so small that the numerical cancellation between the two contributions 
becomes unstable. Fig. 5 shows that the 0{a s ) W + jet or £ ± + ^™ lssmg _|_ o jet cross 
section is essentially independent of s m i n over a wide range of s m i n . In general, one wants 
to choose the largest s m i n possible, to minimize the running time of the program. In this 
case, a reasonable value to choose is s min = 10 GeV 2 . (One must be careful to note that for 
certain distributions, in particular infrared-sensitive ones, a smaller s m ; n might be required 
for some values of the relevant kinematic variables.) 

One quantity of interest is the dependence of the W + jet cross section on the choice of 
experimental cuts. In principle, including higher orders mimics more accurately the correct 
dependence. At leading order, with 'standard' CDF cuts, 

a(W + jets) = 0.781°;°* nb, (4.19) 

This does not depend on the jet defining cut, E^ min , si nce at leading order there is no parton 
in the final state. At next-to-leading order, this is no longer true and in Fig. 6 we show the 
next-to-leading order W + jet cross section as a function of E^ in for the same range of 
scales. As E^ min becomes large, this cross-section approaches the inclusive W cross section. 



B. H X H 2 -» V + 1 jet 

We now turn to vector boson production in association with a single jet. As before, 
the structure of the next-to-leading cross section is described by Eq. |2.12| , however the 
contributing parton-level processes include both those of Eq. |4.12| as well as, 

qq -> V + qq } qq -> V + gg, qg -> V + qg, gq -> V + gq, gg -> V + qq. (4.20) 



As before, the lowest order matrix elements for Eq. |4.20j are obtained by crossing the matrix 
elements given in Appendix A of [15]. 



The finite next-to-leading order matrix elements are given by, 
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(4.21; 



where /C(Qi; 1; Q2) * s given by Eq. [4 . 9| with n = 1. The finite one loop contributions are given 

2 

by Eqs. A.42-A.46 of [|l^]. Once again, it is trivial to cross both /C and «S M (Qi; 1; Q2W 1 ■ 
However, some care must be taken in crossing T since although crossing the helicity structure 
in Eq. A. 43 is straightforward, the coefficients a;, $ and <5j in Eq. A. 44 are expressed in terms 
of the function R(x,y) (Eq. A. 45) which has been written assuming that < x, y < 1. For 
crossed processes this is no longer true. For example, when, 



x < 0, y < 0, 



(4.22) 



in which case, 



R(x,y) 



Alternatively, if, 



-Li 2 



. 1 2 7T 2 

x) -^\og{l-y) + — 



X 



Li 2 



1 



(4.23) 



x<0, y>l, 



(4.24) 



in which case, 



1 X 

R{x, y) = -- log 2 (i - x) + - log 2 (y) + log 



— Li 9 



1 



1 — x 



— J log(y) + y 



Li 2 



(4.25) 



We have checked that crossed finite virtual contributions agree with the results of Gonsalves 
et al. when the vector boson decay current is replaced by its polarisation vector. 

With these crossed matrix elements, we can construct a Monte Carlo program to numer- 
ically evaluate the fully differential vector boson plus one jet cross section at next-to-leading 
order. As before, the vector boson decays are easily included which allows experimental cuts 
to be placed directly on the observed leptons. To demonstrate that the W + 1 jet cross 
section is essentially independent of the unphysical parameter s m j n , Fig. 7 shows that the 
0(a 2 s ) W + 1 jet cross section as a function of s m i n . For s m i n in the range 1-10 GeV 2 , the cross 
section is not dominated by systematic errors and does not depend on s min . We therefore 
set s min = 10 GeV 2 as in the W + jet case. 



An important issue in W + jets events is the signficance of corrections to leading-order 
results. For most quantities, radiative effects should be small so that we can rely on leading 
order to describe the basic features of the data even though the overall normalization of the 
cross section is uncertain. One of the most fundamental distributions is the jet transverse 
momentum distribution which is shown in Fig. 8. 

For the leading order results we choose two renormalisation scales, = My// 2 and the 
total invariant mass of the event, [Ir = s/2. The first scale is the smallest scale available 
and generates the hardest transverse momentum distribution while the second scale is the 
hardest scale available and leads to the softest transverse momentum spectrum. The band 
defined by these two scales represents the range of leading order predictions. The factor 
1/2 is present in the choice of scale so that the corresponding total cross sections at leading 
order (0.117 nb and 0.106 nb) are close to, and bracket, the next-to-leading order result of 
0.113 nb, which is essentially independent of the renormalization scale (here taken to be 
M w ). As can be seen in Fig. 8, the next-to- leading order distribution is somewhat softer 
than the leading-order results. In fact it is even softer than leading order with the largest 
scale. This implies that the standard jet algorithm does not take into account an important 
effect thereby leading to large radiative effects for high transverse momentum jets. 

From a physical point of view, it is clear why the spectrum softens more dramatically 
than can be expected from leading order with the standard jet algorithm. For high transverse 
energy jets, the accompanying soft radiation increases with the energy of the jet. Therefore, 
with a fixed transverse momentum cut ( in this case -Ej^in = 15 GeV), it is easier for the 
soft radiation in the event to fluctuate so that it passes the minimum transverse momentum 
threshold and subsequently be counted as a extra jet. Since we are looking at the exclusive 
jet cross section this event will be removed from the W + 1 jet cross section and added to 
the W + 2 jet cross section. This effect gets more severe when the jet transverse energy gets 
larger leading to a depletion of the W+l jet cross section and a softening of the transverse 
momentum distribution. At leading order this effect is not modelled at all because the 
leading order prediction associates all the energy with the jet and allows no soft radiation 
outside the jet cone. In contrast, at next-to-leading order the hadronic energy around the 
jet cone is modelled, allowing the generation of softer jets within the tail of high transverse 
energy jets. 

From a more mathematical point of view, it is also clear what happens in the exclusive 
jet cross section with high transverse momentum jets. From Eq. (4.9), we can see that 
the high transverse momentum jets generate correction terms of order — a s \og(E } ^ t / E 3 ^ min ) 
which become large if the is much larger than the minimal transverse energy. This is 
undesirable since it implies large radiative effects which are due entirely to the jet algorithm 
itself. In principle, the jet algorithm should minimize these effects in order to be able to 
compare theory with experiment. 

This requires a slight modification in the jet algorithm. By scaling the minimal trans- 
verse momentum cut with the hardness of the event (e.g. summed scalar energies or total 
invariant mass), we allow the hard jets to radiate accompanying soft energy without gener- 
ating additional small jets. That is, we allow the jet to 'vent' its energy without producing 
a large number of soft jets. Events which formerly contained W + 2 jets where one of the 
jets is relatively soft are now counted as If + 1 jet events, thus increasing the W+l jet 
cross section at high transverse momentum. For example, by demanding E^ in = \y/§, the 



correction term is — a s log 2 (A) where we can now choose the constant A and thus control the 
size of the corrections. 

To demonstrate this effect we show in Fig. 9 the factor which we need to multiply the 
leading order distribution (with /i R = s/2) to obtain the next-to-leading order result. For 
the fixed E 3 ^ in , we see a sizeable correction that depends strongly on the jet transverse 
energy. However, if we take a scaling E^ min = max(15 GeV, 0.1 y/§), we get a result very 
close to leading order with only a small enhancement for soft transverse energy jets. For 
the high transverse energy jets the next-to-leading order prediction is well described by the 
leading order result. 



FIG. 5. The s m j n dependence of the W + jet cross section for 'standard' CDF cuts; 
> 15 GeV, E l T > 20 GeV, E™ iss ' mg > 20 GeV, \r]i et \ < 2, \rj £ \ < 1 and a jet cone size 
AR = \J Acp 2 + At] 2 < 0.7. The structure functions are set Bl of |34|] while the factorisation and 
renormalisation scales are fMp = /j,r = My/- For input parameters we choose M\y = 80 GeV, 
T w = 2 GeV, sin 2 9 W = 0.23 and a s (M w ) = 0.1108. 



FIG. 6. The NLO W + jet cross section as a function of the jet defining cut E^ nin 
= hr = 2M W , M w and M w /2. 



FIG. 7. The s m i n dependence of the W + 1 jet cross section for 'standard' CDF cuts. The 
structure functions are set Bl of |}4j while the factorisation and renormalisation scales are 
= VR = M w . 



FIG. 8. The 0(a s ) and 0{a 2 s ) jet and W boson transverse momentum distribution for 'stan- 
dard' CDF cuts. 



FIG. 9. The ratio of next-to-leading order to leading order jet transverse momentum distribu- 
tions for a fixed E^ min = 15 GeV cut and for a scaled E^ min = m a x (15 GeV, 0.1 x cut as a 
function of the jet transverse momentum. 



V. CONCLUSIONS 



The main theme in this paper is the extension of the general method of ref . |L5| , for dealing 
with final state collinear and infrared divergences, to include partons in the initial state. For 
final state partons, the soft and collinear divergences from the bremstrahlung process are 
isolated using a parton resolution parameter s min . These divergences are proportional to 
lowest order matrix elements and can be combined directly with the divergences from the 
virtual graphs to give a finite cross section, Eq. p.4j , which depends on a dynamical /C factor 
multiplying the lowest-order term, along with a finite one loop contribution T . In order 
to extend this to incorporate initial state partons, we have extended the tree level concept 
of crossing to next-to-leading order processes. This is achieved by (a) through the analytic 
continuation of the dynamical factor K, (Eq. |4.9| ) and the finite one loop contribution T 
(see Sec. [TV] ) into the physical region and (b) through the introduction of universal crossing 
functions, Eq. [2.13| , which are essentially convolutions of the structure functions with the 
Altarelli-Parisi splitting functions. Together, /C, T and C(x) form a set of finite building 
blocks with which one can calculate next-to-leading order cross sections. This is summarized 



in Eq. |2.12| ; because the ordered factorization of the soft poles and the factorization of the 
collinear poles, is independent of the hard process, this equation represents the general cross 
section for any hadronic process. 

As an explicit example of this method, in Sec. [IV], we have taken the next-to-leading 



order matrix elements relevant for e + e — > 2, 3 and 4 partons QT5J and crossed two of the 
partons into the initial state to obtain the cross section for, 

pp -> W ± /Z + 0, 1 jets -> ££ + 0, 1 jets, (5.1) 

at next-to-leading order. The phase space is evaluated numerically with the constraint that 
all \sij\ > s m i n and all final state lepton correlations are retained. This makes it possible 
to implement jet algorithms, detector acceptance effects, and other constraints numerically, 
yielding a very flexible Monte Carlo programs as we discussed in Sec. fVJ. One should verify 
that the cross section is independent of the unphysical parameter s min . For our Monte Carlo 
simulations this is indeed the case, see Figs. 5 and 7. 

It is important to note that throughout this paper we have discussed cross sections that 
are exclusive in the number of jets. As a result, we are interested in calculating the W + 
jet cross section, rather than the inclusive W cross section which can only be identified with 
the W + jet cross section as the mimimal transverse energy cut of the jet becomes very 
large. Similarly, we study the W + 1 jet cross section rather than the transverse momentum 
distribution of the W boson which is not directly measurable. As shown in Fig. 8, this is not 
the same as the jet transverse momentum distribution at next-to-leading order. At large Ej! , 
the next-to-leading order jet pt distribution is significantly softened. This is because the 
existing jet algorithm generates an artificially high jet multiplicity in events containing a very 
hard jet by restricting the hadronic radiation around the primary jet in the exclusive jet cross 



section. By modifying the jet algorithm as described in Sec. [TV], these large radiative effects 
can be removed. A detailed study of the implications of next-to-leading order corrections 
to vector boson production in association with 0, 1 jets at Fermilab energies is currently in 
progress p6|j37|] . 

The method we have presented here considerably simplifies the structure of next-to- 
leading order QCD corrections to hadronic processes. It also makes comparison with ex- 
periment more direct through the use of Monte Carlo simulations. Once technical problems 
associated with five point loop diagrams are solved, it should be posssible to use these meth- 
ods to compute processes such as e + e~ — > 4 jets, its crossing pp — ► W ± /Z + 2 jets and 
pp — > 3 jets at next-to-leading order possible. These multijet cross sections are important 
for experiments at LEP and Fermilab because event rates are high and can be studied in 
great detail. 
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